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1 '^ ' Abstract 

f^ ' We examine the regularized zero-range model in an application to 

three-fermion systems - the triton and the hypertriton. We consider 
bound states and low-energy neutron-deuteron and lambda-deuteron scat- 
tering. The model is shown to provide an adequate quantitative descrip- 
^ , tion of these system on a par with finite-range potential models. The well 

r^ . known correlation between the doublet nd scattering length and the tri- 

ton binding energy (Phillips line) finds a natural explanation within the 
model. 
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^ ; 1 Introduction 

The zero-range potential^ is a useful concept which employs separation of scales 
. in physical systems |l|. The model allows selection of the relevant degrees of 

1^^ ' freedom and provides qualitative and often quantitative descriptions of intricate 

S . physical phenomena in an intuitive and transparent way. Conceptually it is 

equivalent to a coordinate space formulation of effective field theories [^ . 

Particularly, the zero-range model effectively reduces a three-body problem 
down to a one-body problem and thus provides a useful means of model inde- 
pendent insights into the important properties of various systems ||3|, y, |5|, |6| . 

However, a direct application of zero-range potentials to three-body systems 
presents a problem - a collapse of the system known as the Thomas effect Q . 
To have a predictive power the zero-range model needs to be regularized. 

Several momentum-space regularization approaches have been recently re- 
ported, which basically amount to a cut-off at high momentum [pi 0|, or a 
combination of a cut-off and an additional three-body force ||l^ . 

A combination of an energy cut-off and an effective three-body force had 
also been used in a coordinate-space problem with ^-interactions pT| . 

Recently we have introduced an alternative coordinate-space regularization 
approach [n2|. Unlike |0| we do not explicitly use (5- functions. The zero-range 
potential is formulated as a boundary condition for the wave-function at the 

^also referred to as "pseudo-potential" or "<5-interaction" 



origin. Regularization is achieved by a suitable generalization of the boundary 
condition and can be thought of as a soft cut-off. 

In contrast to momentum-space theories P, 0, Effl our approach does not 
require an additional three-body parameter to set the scale of the three-body 
binding energy. 

The purpose of this paper is to test the suggested regularized zero-range 
model in an application to systems of fermions, considering both bound state 
problems and scattering processes, and comparing, where possible, with finite- 
range potential models and effective field theories. 

We shall give a brief introduction to the regularized zero-range potential 
model, consider the hyper-spheric formalism of "1-1-2" scattering, and finally 
apply the model to the nnp and Anp systems. 

In the following we shall refer to the regularized zero-range potential as 
contact interaction. 



2 Three-body problem with contact interactions 

In this chapter we discuss hyper-spheric adiabatic solution of the three-body 
problem with contact interactions. We use the definition of the two-body 
Schrodinger equation with contact interaction as a free s-wave equation with 
the solution ipk{r) = shi(kr + S)/r, and the boundary condition at the origin 

-^^rijjk ^kcot5^- + -Rk^ + PRk'^ , (1) 

ripk dr ^^g a 2 

where 5 is the two-body phase shift, a, R, and P are the scattering lengthR 
the effective range and the shape parameter and k is the wave- number. The }? 
and fc"* terms are necessary for the contact interaction model to have a regular 
three-body ground state solution Q . 

2.1 Hyper-spheric adiabatic wave- function 

We use the hyper- spheric adiabatic approximation [|l3| which employs the hyper- 
spheric coordinates p and 17 (see the appendix). The total wave- function ^P of a 
three-body system is written as a product of a hyper- angular function ^(p, il) 
and a hyper-radial function f{p) 

^^p-"^I{p)HP.^)- (2) 

The function $(p, Vt) is an adiabatic eigenfunction of the three-body Hamilto- 
nian with a fixed hyper-radius p. The corresponding eigenvalue \{p) serves as 
an effective potential for the hyper-radial function /(p) 

where E is the total energy of the three-body system, m is the mass unit used 
in the definition of p (see the appendix) , and 



/o2 
d^'^*{p,n)—<^{p,^). 



(4) 



•^notc the sign convention for the scattering length in 



In most cases the one-channel adiabatic approximation (0) proves to be quite 
accurate WM. However, for higher accuracy several adiabatic eigenf unctions 
might need to be included. Additional cigcnfunctions are also needed, for ex- 
ample, for a description of break-up processes. Although inclusion of additional 
eigenfunctions is straightforward, in the present context we restrict ourselves to 
the lowest adiabatic channel only. 

2.2 Adiabatic eigenvalue equation. 

We assume that contact interactions act only on the s-waves and vanish every- 
where except at the origin. The adiabatic eigenfunction <&(p, fi) can be then 
suitably expressed in terms of the free s-wave Faddeev components ip^{ai) 
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(6) 



where v = VA -I- 4, and at is the i-th hyper-angle. In the following with shall 



refer to both v and A : 



4 as angular eigenvalues. 



The quantities Ai{p) and i/(p) are to be determined from the boundary 
conditions (11]). The latter can be conveniently reformulated in terms of the 
hyper-angles as |l2| 
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where a^, Ri, Pi and fXt are the scattering length, effective range, shape parame- 
ter and the reduced mass (divided by the mass unit m) for the two-body system 
of particles j and k. We assume here and in the following that the indices 
{i,j, k} form a cyclic permutation of {1,2,3}. 

For small ai the wave-function $ can be expanded as ||l3|] 



sin(2ai)<l> 



Ai{p)(p^{ai) 



+ Aj{p)2a., 



sin(20ji) 



Ak{p)2a, 



sin(20fej) 



0{a1 



(f>ji = arctan 

The boundary conditions 
for the coefficients Ai{p) 



lmk{mi +m2 + ma) 



(8) 
(9) 



mjirii I 

then turn into a set of three linear equations 
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E 



M,jAj{p) =0, 

M.,=2^^,.^, 
sm(20jj) 



(10) 



where </3'i,(0) = dLpi,{a)/da\^^i^ and where we have used a short-hand notation 
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The set of equations ( pJ]| ) has nontrivial solutions only when the determinant 
of the matrix M{p, v) is zero, which gives the equation to determine i'{p) and 

detM(p,i/) =0. (12) 

3 "1+2" scattering 

In this chapter we discuss the hyper-spheric formalism for the elastic scattering 
of a particle against a bound system of the other two. This process is often 
referred to as "1-1-2" scattering. 



3.1 Adiabatic solution at large distances 

A bound two-body subsystem - say, number 3 - with the binding energy B — 
h K^/ {2mix^) gives rise to an eigenvalue v{p) which at large p asymptotically 
behaves as lll2l 



. P 



/P-3 

2mB 



Kp) 
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n 
with the corresponding (normalized) adiabatic wave-function (|6|) being 

PK3 



^v{a:i) == ./2^=^exp 



as 



(13) 
(14) 

(15) 



Indeed, (p^{(t>2i) vanishes exponentially for large p and the 3-rd equation in 
(R) asymptotically decouples from the other two and turns into 



as 2 



/psi^ 



p^m 



/P-iV 



(16) 



which with the substitution v 
for the two-body bound state 



P J \ P 

^pI \fpi takes the form of the usual equation 



ik — k cot 6{k) 



(17) 



with the solution iK^. Note that since ^Pu{4''ii) vanishes exponentially with p the 
correction to (03) is also exponentially small p2[ . 

The additional term Q in the hyper-radial equation (ph asymptotically ap- 
proaches Q -^ — l/(4p^), which can be directly verified from dla). Although 
this term is generally small it is important to ensure the correct asymptotics of 
the hyper-radial equation in the case of a bound two-body subsystem |l4 . For 
simplicity we shall in the following always use only the leading term — l/(4p^). 



often referred to as the Langer correction term. The hyper-radial equation then 
becomes 

with the correct "1+2" asymptotics at large p 

d^ 2m 



2 iE + B)]fip) = 0. (19) 



3.2 Scattering wave- function and phase shift 

Within the contact interaction model the two-body bound state wave-function 
simplifies to exp(— K3ri2)/ri2, where ri2 is the relative coordinate. The wave- 
function describing an elastic scattering of the third particle against the bound 
system of the first two must then have the following asymptotic form 

* ^ — exp(-K3ri2) sin(fcr(i2)3 -l- 6) , (20) 

ri2 '^(12)3 

where ?'(i2)3 is the relative distance between the third particle and the center 
of mass of the first two, 6 is the sought "1-1-2" phase shift, and the scattering 
energy is ?i^fc^/(2m^3). 

We shall show now that the three-body wave-function (0) indeed takes this 
asymptotic form at large p. The asymptotic form ( |l9[ ) of the hyper-radial equa- 
tion has a solution 

fip)^sm{qp + A), (21) 

where A{q) is the hyper-radial phase and g^ = 2m{E + B)/h > 0. Multiplying 
the hyper-radial function (p^) and the adiabatic function (RS) gives the total 
wave-function 

vP = p-^/2sin(gp + A(g))— L^ /^expf-^a3) . (22) 

sm(2a3)Y V^ V VMs / 



For pn^/y/JTi 3> 1 the exponent in (22) is non- vanishing only when as <C 1 



Therefore in this region pas — x and p — y (neglecting terms of the order of 
a§) and the wave- function (p2h becomes 



*~-exp( '^x] - sm{qy + A{q)) . (23) 

^ \ VJJ^ J y 

Recalling that x — ^JJiirn and y = ^//^(i2)3''(i2)3 (see the appendix) we obtain 
^ ~ — exp (-K3ri2) sin(fcr(i2)3 + A(g)) , (24) 

■ri2 7-(i2)3 

where fc — g^M(i2)3- This is exactly the sought asymptotic form (|2^). The 
hyper-radial phase shift Ai^q) is then equal to the "1-1-2" phase shift b(k). 

Thus the " 1-1-2" phase shift ^(fc) can be easily calculated within the hyper- 
spheric approach by solving the hyper-radial equation (^ with the wave-vector 
9 = fc/^//i(i2)3 and extracting the hyper-radial phase A(g) from the asymptotic 
form (pl|). The "1-1-2" phase shift ^(fc) is equal to the hyper-radial phase A((7). 
The " 1-1-2" scattering length can then be found as the limit a^^ — limfe^o ^ cot 8. 



Table 1: Parameters of the contact interaction (m in the ^Si (triplet) np and 
^5*0 (singlet) nn channels. Here a is the scattering length, R - the effective range, 
P - the shape parameter, Ed is the energy of the bound np state (deuteron), 
and Eyirt is the energy of the virtual (anti-bound) nn state. 



VnpCSi) 


a, fm 


i?, fm 


P 


Ed, MeV 


I 

II 
III 


-5.26 
-5.14 
-5.02 


1.54 
1.38 
1.21 


0.084 
0.143 
0.256 


-2.2 
-2.2 
-2.2 


Kn(^5'o) 


a, fm 


R, fm 


P 


Ey„t, MeV 


I 

II 
III 


19.00 
19.09 
19.21 


2.10 
1.91 
1.64 


0.086 
0.132 
0.238 


-0.1 
-0.1 
-0.1 



4 Contact interactions for np and nn systems 

The basic quantities that characterize the np interaction at low energies are the 
phase shifts and the deuteron binding energy. We therefore choose the contact 
interaction parameters a, R, and P in (|l|) by fitting the phase shifts while 
keeping the deuteron binding energy fixed at -2.2 MeV. For the nn system we 
instead keep the energy of the virtual s-state at -0.1 MeV. Our parametrization 
does not allow change of sign of the phase shifts at large energies. Therefore 
we can only aim at describing the phases below 100 MeV for nn and below 
150 MeV for np scattering. 

For each pair of particles we consider three sets of parameters labelled I, 
II, and III (see Table 1), where II is more or less the best possible fit, while I 
and III represent small variations towards better description of correspondingly 
lower and higher energies. 

The np and nn phase shifts for these models are shown, correspondingly, on 
Fig. |l| and Fig. § We have used ^0=197.3 MeV fm and m^^nip =939 MeV. 
Potentials I provide a better description at lower energies while overestimate the 
phases at larger energies. Potentials III are better in the higher energy region 
at the expense of the lower energies, while potentials II represent a reasonable 
compromise. 

The larger P-parameters give smaller phase shifts at higher energy or, cor- 
respondingly, weaker interaction. One can say that the P-parametcr smoothly 
cuts off the interaction at higher energies, thus regularizing the three-body prob- 
lem in. 



5 Application to the nnp system 

The spins of three nucleons can couple to a total spin S of 3/2 and 1/2, which for 
relative s-waves is also the total angular momentum. In the S =1/2 (doublet) 
case there is a bound state of a proton and two neutrons - a triton. No three- 
body bound state is found in the S—3/2 (quartet) case. We assume that the 
neutron and the proton have the same mass ?7i=939 MeV, which is also used as 
the mass scale in the definition of p (see appendix). 
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Figure 1: The '^Si np phase shift S as function of the center-of-mass energy 
Ecm for the contact interactions from Table 1 . The circles represent the results 
from the Argonne vis potential |l3] . 



nn phase shifts 




Er. 



60 80 

(MeV) 



Figure 2: The ^Sq nn phase shift S as function of the center-of-mass energy 
Ecm for the contact interactions from Table 1 . The circles represent the results 
from the Argonne tiig potential Q . 



5.1 Quartet nd scattering 

5.1.1 Adiabatic eigenfunction 

In the case of particles with spins we have to include also the spin-functions C 
in the adiabatic wave-function 



^p,n) = J2QA 



' shi{2ai) 
For S'=3/2 the spin functions C,i are represented by 



(25) 



(26) 



where ti designates the spin-function of the «-th particle where the z-projection 
of the spin operator is equal to +1/2. 

Let the proton be particle number 1 and the neutrons number 2 and 3. Two 
neutrons with total spin 1 can not be in an s-state. Therefore the Faddeev 
component number 1 must vanish identically. Under exchange of the neutrons 
the angular functions fu{'^2) ^nd 99 j, (0^3) in (Ea) transform into each other while 
the spin functions Q do not change. Therefore the antisymmetric wave-function 
for S =3/2 is 

Vu{a2) V^{az) 



$(p, Vt) = A2C2 



sin(2a2) sin(2a3) 



(27) 



5.1.2 Eigenvalue equation 

The system of three equations (lOh reduces in this case to only one equation 
which includes triplet neutron-proton parameters 



where we have used that (ha : 



4 TT 



Lnp^ui^) , 



(28) 



7r/3 for three particles with the same mass and 



where L„p, see eq. (|l l|) , contains triplet neutron-proton parameters a„p, i?„ 
and P, 
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(29) 



Using the explicit expressions (|g) for the functions ipi, the eigenvalue equation 
281) can be written as 

i^cos(j/f) + ;^sin(z/f) 



sin(z^l) 



P 






P R'^ 

^ npJ-^np 



y/JIv 



(30) 



where we have dropped the index at the reduced mass fj, as we assumed that all 
three nucleons have the same mass. 



Quartet nnp adiabatic eigenvalue 
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Figure 3: Quartet adiabatic eigenvalue A as a function of p for the three np 
interactions from Table 1. 



For small p this equation has a root which to lowest orders in p is 



64{^f^PnvR^ 



,3 



A = 12- 



8(V^) 



3 jr-p /?3 



(31) 

(32) 



The solution A(0) = 12 is an ordinary solution for regular finite-range potentials 
[g_3| corresponding to the hyper-angular quantum number K—2. 

When Rnp and Pnp are equal to zero, as in the ordinary zero-range potential, 
the eigenvalue equation ( |30| ) at p <C Unp has a real root f =2.166. Therefore 
there are no Thomas or Efimov effects in the quartet case even without regu- 
larization [Q. 

For large distances, due to the bound two-body subsystem (deuteron) the 
eigenvalue is asymptotically approaching 



. P 
2mBdP^ 



4, 



(33) 
(34) 



where Bd— 2.2 MeV is the deuteron binding energy and Kd 
(Kd =0.23 fm-i). 



2'mpBd/h^ 



5.1.3 Adiabatic eigenvalues and results 

The solutions \{p) of the equation ( pO| ) for the three np potentials from Table 1 
are shown on Fig. 0. The difference between the curves is not big and, moreover, 



Table 2: Quartet nd scattering length, ^Und- Here Vnp refers to the np contact 
interaction from Table 1. 



*a„d, fm 



Vr 



UT 



np 



Vr 



WT 



np 



Vr 



TTTIY 



np 



-7.03 -6.86 



-6.67 



EFT |16| 



-6.33 



Exp. Il'i 



-6.35±0.02 



the region, where the curves differ appreciably, is largely limited to p <5 fm. 
However, in this region A -|- 4 >0 and, consequently, the effective potential 
(A -I- 4)/p^ in the hyper-radial equation (]l8) is repulsive. Therefore the hyper- 
radial wave-function in this region is small and thus insensitive to the small 
variations in the eigenvalues A. 

The quartet nd scattering lengths calculated using these three adiabatic 
eigenvalues in the hyper-radial equation (^8|) are given in Table 2. All three 
contact interaction models give similar results which are very close to the ex- 
periment and to the result of Bedaque and van Kolck Q obtained within the 
effective field theory approach. 



5.2 Doublet nd scattering and the triton 

5.2.1 Adiabatic eigenfunction 

For the doublet (5=1/2) state we have to include in the adiabatic function ( p5|) 
three different spin- functions, (T„„, r^p, Tps where the two neutrons, ni and n2, 
are in the singlet state while the neutron-proton pairs, n2P and pni are in the 
triplet state 



((ni (Xin2)o <X)p)i 



o (TsTl) i2 



2P = (("2<8)p)i <X"T-l)l 



Tp3 = ((p«)ni)^«)n2)i = \/ g (T1T2) is 



V2 



(T2i3 - i2T3) Tl 



^ (Tail + isTi) T2, 



^ (Tli2 + ilT2) T3, 



(35) 



where {a ^ b)s denotes coupling of the spins of the particles (or systems of 
particles) a and b to the total spin s. 

The wave-function must be antisymmetric under exchange of the two iden- 
tical neutrons, which in our case are particles number 2 and 3. Under transfor- 
mation 2< — >3 the wave-functions undergo the following transformations 



^nn ^ ^nn-i '^2p ^ ^ ^p3 1 



(36) 



where ^:ll,{a.i) are the usual solutions (||) to the free Faddeev equations. The 
antisymmetric hyper-angular function can therefore be written as 



smzai Sin2a2 



Tp3^2 



sin 2a3 



(37) 
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5.2.2 Eigenvalue equation 

The hyper-angular function must satisfy the spin-projected boundary conditions 



dai 
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9a2 



a2 



(r.» 



L„ 



Q2=0 



a2 



(rl^) 



Q2=0 



(39) 



where L„„ and i„p (see eq. (|ll|)) contain, respectively, the singlet neutron- 
neutron (a„„, Rnn, Pnn) and triplet neutron-proton (a„p, i?„p, P„p ) parameters 
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(40) 
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1 + t^npti^^ 


V P / _ 


(41) 



Using the explicit form (^) leads to a system of linear equations for the coeffi- 
cients Ai 



^I'/'UO) + ^2 (fTL^2p) 
= Lrm^l<y5i/(0), 
^2¥'U0) + ^1 (^sVnn) 



4^.(^3) 
V3 


^ /.t _ \ 4</'i.(7r/3)' 

- ^2 «„7-p3J 




4^.(7r/3) 
V3 


^2 (tzVps) /^^'^'(t^/S) 



(42) 



(43) 



Again this system of equations has nontrivial solutions only when the determi- 
nant of the corresponding matrix M is zero, which gives the eigenvalue equation 



det M = , 



(44) 



Ap'^nnj ^^^^(f) 



V'M - (t-2Vp3) isM'l) - LnpMO) 



where the cross-products of the spin functions can be easily calculated from (p5|) 



I _ _\/3 I _ v/3 I 



(45) 



Without the regularizing terms, that is when effective ranges and shape param- 
eters are zero, the eigenvalue equation ( [f4|) has in the limit p <C a an imaginary 
root v = 0.787i which results in the Thomas collapse and the Efimov effect ||l8|| . 

With non-zero effective ranges and shape parameters the equation ( [f4| ) has 
a real root h'{p = 0)=0 and the Thomas collapse is therefore removed. The 
three-body system then has a well defined ground state. The Efimov effect can 
still exist in the limit a — > oo. 

The large distance asymptotics is determined by the deuteron binding energy 
and therefore is the same as for the quartet case, eq. (B3h. 
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Doublet adiabatic eigenvalue 




Figure 4: Doublet adiabatic eigenvalue A as function of p for different combi- 
nations of contact interactions from Table 1. 

Note that our regularization approach, contrary to M BJ and flG], utilizes 
only two-body parameters and does not require an additional three-body pa- 
rameter. In principle, we could also introduce a three-body contact interaction 
in the hyper-radial equation and use this additional parameter to further fine 
tune the model. However, this is not the purpose of this paper. 



5.2.3 Adiabatic eigenvalues and results 

Contrary to the quartet eigenvalue, which only depends on the np parameters, 
the doublet eigenvalue depends upon both Vnp and Vnn- Fig. 4 shows the 
solutions of the eigenvalue equation (H4) for all possible combinations of Vnp and 
Vnn from Table 1. The doublet effective potential (A -I- 4)/p^ has an attractive 
pocket where variation is somewhat larger than in the quartet case. One also 
notices that all curves have basically the same shape, the only difference being in 
the depth of the attractive pocket. This property is essential for the discussion 
of the Phillips line in the next subsection. 

Numerical solution of the hyper-radial equation (RSh with these eigenvalues 
gives in each case one bound state, the triton, the energy of which is shown in 
Table 3 together with the corresponding doublet nd scattering length. Only the 
results for the corresponding potentials are shown in the table. Other combina- 
tions give the results which fall somewhere in between the shown results. 

One can conclude that the contact interactions adequately describe the triton 
and the low-energy doublet nd scattering. The potentials II even provide a good 
quantitative description similar to finite-range potential models. 



5.2.4 Phillips line 

The Phillips line ||l^ is an established relationship between the triton binding 
energy and the doublet neutron-deuteron scattering length - different models 
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Table 3: The triton ground state energy, Et, and the doublet nd scattering 
length, ^a„d, for the nucleon-nucleon interactions from Table 1. 



Vnn 


^7ip 


Et, MeV 


^a«d, fm 


I 

II 
III 


I 

II 
III 


-11.4 
-7.25 
-4.98 


2.82 
-0.61 
-3.15 


Exp. 


-8.48 


-0.65 
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Figure 5: The Phillips plot of the doublet nd scattering length ^a„d versus 
the triton binding energy Bt- The points indicate the results from different 
combinations of the contact interactions form Table 1. The solid line represents 
the finite-range potential models |20 . The dashed line represents the results 
from the renormalized effective field theory with three-body force pj . 



for the nucleon-nucleon interactions with the same low-energy properties pro- 
duce different but correlated values for these two quantities. Graphically this 
correlation is customarily represented as a curve (Phillips line) in the plot of 
the quartet nd scattering length '^and versus the triton binding energy Bt- 

The Phillips line for our model is shown on Fig. 5 where we have plotted 
the results for different combinations of the np and nn interactions from Ta- 
ble 1. Also shown are the results from finite-range models ^^ and from the 
renormalized effective field theory with three-body force [El| . 

Our Phillips line is reasonably close to but still somewhat different from 
effective field theory curve. This indicates that the analytic properties of these 
models are different. 

Now the question to answer is why do the results fall on a one-parameter 
curve while we in principle vary at least two independent ingredients - the 
nn and np interactions? In our model the answer can be most easily seen on 
Fig. 4 where shown are the eigenvalues for all combinations of the nn and 
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np interactions from Table 1. Due to the analytic properties of the eigenvalue 
equation (Q) all curves start from A(p = 0) = and end up at A(p » i?) = 
—2mBdP^h — 4. The shape of the curves is largely preserved - there are no 
crossings (with one exception) - and the difference is practically only in the 
depth of the attractive pocket. Therefore we have got basically a one-parameter 
family of the adiabatic potentials, the parameter being, for example, the depth of 
the pocket. The one-parameter family of adiabatic potentials would necessarily 
result in a single curve in the Phillips plot where this parameter is changed 
along the curve. 

There are in fact two curves on Fig. 4 which, despite being almost identi- 
cal, do cross. Correspondingly these curves produce two close points (around 
6.3 MeV) on the Phillips plot. Fig. 5, which are slightly off the line. 

Note that in the effective field theory approach it is the three-body cut-off 
parameter that is changed along the curve ||2l|| . In our language that corre- 
sponds to fixing the two-body parameters, thus fixing the adiabatic potential, 
and instead varying the strength of the three-body force (one parameter). This 
would, of course, also lead to a single curve on the Phillips plot. 

6 Application to Anp system 

The Anp system has a weakly bound ground state, hypertriton, with the spin 
1/2 and the binding energy of 0.13±0.05 MeV relative to the deuteron. Many in- 
vestigations have been carried out with different interaction models (see ||2^, ^ 
and references therein) some of them including A-E conversion |23] . The results 
exhibited considerable variations mostly due to uncertainties in the lambda- 
nucleon interaction. 

One might expect that a certain correlation exists between the Ad scattering 
length and the hypertriton binding energy, similar to the Phillips line for the 
nnp system. The purpose of this chapter is to establish the Ad Phillips line for 
our contact interaction model and compare it with the finite-range model. We 
shall ignore the A-E conversion which is beyond the scope of this paper. The 
mass of the A particle is assumed to be 77ia=1116 McV. 

6.1 Contact interactions 

For the np subsystem we shall use the interaction II from Table 1. Due to 
lack of experimental scattering data we cannot use the same procedure for the 
AA^ subsystem. Instead we simply fix the scattering length and the effective 
range and vary the shape parameter around the values similar to those of the 
nn interaction. In |g2| a set of low-energy parameters was found which gives 
a reasonable description of the hypertriton within the finite-range three-body 
model. We assume that pA and nA interactions are identical and adopt the 
scattering lengths and effective ranges in Table ^. They are close to the set Elr 
in |22[| . The singlet Pg and triplet Pt parameters are for simplicity assumed to 
be equal and are varied in the range 0.08-0.11. 
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Table 4: Parameters of the ^5*0 and '^Si A-nucleon contact interactions. The 
singlet and triplet shape parameters are assumed to be equal and are varied in 
the range 0.08-0.11. 





a, fm 


R, fm 


P 


'So 


1.9 


3.5 


0.08-0.11 


's. 


1.5 


3.7 


0.08-0.11 



6.2 Adiabatic eigenf unction and eigenvalue equation 

We assume that the np pair is in the triplet state while pA and nA pairs can 
be both in triplet and singlet states. We shall therefore use the following spin- 
functions 



' np 
CpA 
TpA 
CriA 

TnA 



((n®p)i®A)i = y^(UTp)u-y^(Trap + inTp)TA, (46) 



(47) 



((p«)A)o(g)n)i = \/ 2 (TpiA - IpTa) Tn, 



Up ® A)i ® n) 1 = ^^ (TpTa) i« -y ^ (TpiA + IpTa) Tn, (48) 



(49) 



((A(8)n)o(X)p)i = -\/-(TAin - iAT«) Tp, 



((A(K)n)i(g)p)j 



^ (TaU) ip -Jl iUln + iAtn) Tp • (50) 



We assume that the wave-function is symmetric under exchange of the two 
nucleons. Since under exchange n *-^ p the spin functions transform as 



fpA *~^ "CTnAj TpA ^^ TnAi Tnp ■*"*■ T„p, 

the symmetric hyper-angular function (E3) can be written as 



(51) 



sm 2a 1 



^''^"^^ {Ba,A + Ct,a) + ^^ {-Bar^A + Ct^a) • (52) 



sin 2a2 



sin 2a3 



For this system we have three spin-projected boundary conditions (R), the triplet 
np, and singlet and triplet A'^A 



d 



^- l^a^'"^'^ 



a, =0 



4*^,lp(2ai'i')a,=0' 
^iA(^2a2<J>)^^^ = 4V^a(2«2$)„^.o^ 



^.'a (a!^^"^'^ 



4*^4 (2«2$), 



where s and t refer to singlet and triplet states and where, as usual 

2 






f^ab \ a 



1 



1 d(c) (v^/JI^ 

2 "^^ l^^ 



P' 



)(c) /p(c 
ab Y^ab 



)\3 /v^/JI^ 



P 



(53) 
(54) 
(55) 

(56) 
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where fiab — {^ / 'ni)ma'mh / {ma + rrii,) is the reduced mass of particles a and 6 in 
units of m. The resuhing system of hnear equations for the coefficients A, -B, 
and C is 



V.(0) + 2^^#i^ (Brtp^pA + CrtpT.A) = ^il>.(0), 



sin 2(/i2i 



sin 20; 



'21 



pA'np 



^ , ,r,x , ^ 2(p^(i^2i) t 



sin 20: 



21 



pA'np 



2y^(023) / p t 



y^ t 



A + CTp^TnA 



C4*lr^^(0). 



And, finally, the eigenvalue equation 



det M = 



(57) 

(58) 

(59) 
(60) 



f ^'M-^mL% 


r,2¥>.(02l) t ^ . 


22v.(02l) t r„. \ 
^ sin202l 'np'pA 


2v„W2l) t , 
sin 2021 >A^"P 


^Uo)-^.(o)4'j^ 

2v'..(te) t _ ^ 

sin 2023 "^pA^nA 


2^.(023) ^t _ ^ 
sin 2023 %A^»A 


2v,.(02l)_t _ 
I sin2</,2l 'pA'"P 


2¥'.(te)_t _ ^ 
sin 2023 'pA^nA 


1 2>p,.(023),t , , 

^" Sin 2023 >A^nA / 



where the cross products of the spin-functions are 



T^pCTpA = 


\/3 ^ 
2 ' ''"p'^pA ^ 


1 t 


V3 
2 


cr^A'^"A = 


1 t 


n/3 
2 ' 




TpACTnA = 


V3 t 
2 ' '^pA'''"A — 


1 
^2' 





(61) 
(62) 
(63) 



6.3 Eigenvalues and results 

The solutions of the eigenvalue equation (pw with parameters from Table 
are shown on Fig. ^. Again the asymptotics at both p=0 and p ^ R are 
independent on the shape parameter. The agreement between the curves in the 
asymptotic region is better than in the triton case because here we keep the 
two-body scattering lengths and effective ranges fixed (Table ||) while in the 
former case they were slightly varied (Table 1). The depth of the attractive 
pocket, however, does show some variation with the shape parameter. 

The corresponding hypertriton binding energies and doublet Ad scattering 
lengths are shown on the Phillips plot. Fig. H, together with the results from 
finite-range models |22]. The agreement is even better than in the triton case, 
although the curve for the contact interactions is still marginally below. Basi- 
cally the contact interaction model proves to be on a par with the finite-range 
model. 

Note that this is not a simple repetition of the triton results as the hypertri- 
ton has a very different structure - the A particle is not identical to the neutron, 



16 



Knp adiabatic eigenvalue 



-< 




' I ' ' ' I ' ' ' I 
P,=Pt=0.08 — 

P,=Pt=0.09 -- 

P,=Pj=0.10 --- 

P,=Pi=0.11 



_L 



1 1.5 2 
p, (fm) 



2.5 



Figure 6: Doublet adiabatic eigenvalue A as function of p for Knp system with 
contact interactions from Table ^. 
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Figure 7: The Phillips plot of the doublet Ad scattering length "^a^d versus the 
triton binding energy relative to deuteron (Bd=2.2 MeV). The curve is obtained 
from the contact interaction model. The circles represent finite-range models 
[0. 
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which results in quite different angular functions (see ( p7\ ) and (p2h), and also 
does not form a bound state with the proton. 

The better agreement reflects also the fundamental property of weakly bound 
systems - halos (M - of which the hypertriton is a remarkable example, to be 
insensitive to the properties of the underlying interaction models. 

7 Conclusion 

We have applied the recently formulated regularized zero-range model to fermionic 
three-body systems nnp and Anp. The model employs low-energy two-body pa- 
rameters, scattering length and effective range, and the shape parameter. The 
three-body problem is solved within the hyper-spheric adiabatic approach where 
the zero-range model allows analytic solutions for the adiabatic hyper-angular 
functions. 

The large-distance asymptotics of the adiabatic potential is completely de- 
termined by the scattering length and effective range and therefore is model- 
independent. The shape parameter mostly determines the depth of the attrac- 
tive pocket and is thus important for the actual value of the binding energy. 

Our regularization approach does not require an additional three-body pa- 
rameter in contrast to the alternative momentum space regularizations. The 
finite triton binding energy is obtained using only two-body data, much like 
with the finite-range models. 

Once the two-body parameters are fit to the nucleon-nucleon phase shifts 
the model provides an adequate quantitative description of the nnp system - 
the bound state, triton, and the low-energy nd scattering. The established 
correlation between the doublet nd scattering length and the triton binding 
energy, the Phillips line, finds a natural explanation within the model. 

For the Anp system, with its weakly bound state, hypertriton, the Phillips 
lines for zero-range and finite-range models are almost identical, refiecting the 
fact that the topical weakly bound systems, halos, are residing mostly at large 
distances, where the adiabatic potential is independent upon the underlying 
potential model. The regularized zero-range model is therefore ideally suited 
for such systems. 

The simple analytic functions of the model might serve as a convenient 
Sturmian basis in other three-body applications. 

In conclusion the regularized zero-range potential model can be effectively 
used for quantitative descriptions of various three- and many-body systems in 
discrete and continuum spectra. 
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9 Hyper-spheric coordinates 

If rrii and r^ refer to the i-th particle then the hyper-radius p and the hyper- 



angles ai are defined in terms of the Jacobi coordinates x^ and y^ as [25 



rrij + TOfe J 



/JmItj ~ rfe) , y,; = ^/JJiJ^ir, 



1 mjiTik 1 miimj+nik) ,„.. 

Mi = , M(ifcV = -^ — [ (64) 

m rrij + mk m nii + rrij + nik 

psin{ai) = Xi , pcos{ai)^yi, 

where {i,j,k} is a cyclic permutation of {1,2,3} and m is an arbitrary mass. 
The set of angles Hi consists of the hyper-angle ai and the four angles Xi/|xi| 
and yi/|yi|. The kinetic energy operator T is then given as 

A = — -i— ^sin(2a0-4-' ^^^^ ^ ^^^ 



sin(2ai) da^ * sin^(ai) cos2(ai) 

where Ix^ and ly. are the angular momentum operators related to Xj and y^. 
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